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Abstract 

The Drazin inverse solutions of the matrix equations AX = B, XA — B and 
AXB — D are considered in this paper. We use both the determinantal repre- 
sentations of the Drazin inverse obtained earlier by the author and in the paper. 
We get analogs of the Cramer rule for the Drazin inverse solutions of these ma- 
trix equations and using their for determinantal representations of solutions of 
some differential matrix equations, X' + AX = B and X' + XA = B, where 
the matrix A is singular. 
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1. Introduction 

In this paper we shall adopt the following notation. Let C"*^" be the set of 
m by n matrices with complex entries and \m be the identity matrix of order m. 
Denote by a and a^. the jth column and the ith row of A G C™^", respectively. 
Let A,j (b) denote the matrix obtained from A by replacing its jth column with 
the vector b, and by A^. (b) denote the matrix obtained from A by replacing 
its ith row with b. 

Let a := {ai, . . . , a^} C {1, . . . , m} and P := . . . , C {1, . . . , n} be 



A" 



denotes the minor of A 



subsets of the order 1 < k < min {m,n}. Then 

determined by the rows indexed by a and the columns indexed by /3. Clearly, 
I A" I be a principal minor determined by the rows and columns indexed by a. 

For 1 < A: < n, denote by 

Lk,n := {a : a = (ai, . . . ,afe) , 1 < ai < ■ • • < "fc < n} 

the collection of strictly increasing sequences of k integers chosen from the set 
{1, . . . , n}. For fixed i E a and j € /?, let 

Ik,m{i} := {a : a e Lk,m,i e a}, Jk,n{.j} := : /3 e £fe,n, j G /3}- 



Preprint submitted to Applied Mathematics and Computation 



January 29, 2013 



Matrix equation is one of the important study fields of linear algebra. Linear 
matrix equations, such as 

AX = B, (1) 
XA = B, (2) 

and 

AXB = D, (3) 

play an important role in linear system theory therefore a large number of 
papers have presented several methods for investigating these matrix equations 
(for example, see [ij -[!]). 

In some situations, however, people pay more attention to the Drazin inverse 
solutions of singular linear systems and matrix equations Moreover, Xu 

Zhao-liang and Wang Guo-rong in proved that the Drazin inverse solutions 
of the matrix equations ((U, ^ and @ with some restricts are their unique 
solutions. The Cramer rule for the Drazin inverse solution of the restricted 



system of linear equations are used in [15|- [17|. The Cramer rules for solutions 
of the restricted matrix equations ([!]), (pi) and (|3]), in particular for the Drazin 
inverse solution, are established in (18i]-|2lj. 

In this paper, we obtain explicit formulas for determinantal representations 
of the Drazin inverse solutions of the matrix equations ([Ij , ^ and ([3]) and using 
their for determinantal representations of solutions of some differential matrix 
equations. The paper is based on principles used in [2^, where we obtained 
analogs of the Cramer rule for the minimum norm least squares solutions of the 
matrix equations (1), (2) and (3). Liu et al. in [^l deduce the new determi- 

(2) 

nantal representations of A^g and the Cramer rule for the restricted matrix 
equation ^ based on these principles as well. Since the Drazin inverse and the 
group inverse A are outer inverses Aj, ^ for some specific choice of T and S, then 
the results obtained in p4| generalize in some ways some results of the paper. 
But we get the more detailed representation of the Drazin inverse solutions, 
and therefore we can used their for determinantal representations of solutions 
of some differential matrix equations. 

The paper is organized as follows. We start with some basic concepts and 
results about the Drazin inverse in Section 2. We use the determinantal repre- 
sentation of the Drazin inverse obtained in and also another determinantal 
representation is obtained in this section. In Section 3, we derive explicit for- 
mulas for determinantal representations of the Drazin inverse solutions for the 
matrix equations (H]), ^ and ([3]). These formulas generalize the well-known 
Cramer rule. In Section 4, we demonstrate their using for determinantal repre- 
sentations of solutions of some differential matrix equations, X' -I- AX — B and 
X' + XA = B, where the matrix A is singular. In Section 5, we show numerical 
examples to illustrate the main results. 
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2. Determinantal representations of the Drazin inverse 

For any matrix A G (^"xn ^[^^i IndA =^ where a positive integer k =: 

IndA = min |rankA'^'+^ = rankA*^}, the Drazin inverse is the unique ma- 
fee7Vu{o} ^ ^ 

trix X that satisfies the following three properties 

1) A'^+iX = A'^; 

2) XAX = X; (4) 

3) AX = XA. 

It is denoted by X = A^. 

In particular, when IndA = 1, then the matrix X in is called the group 
inverse and is denoted by X = A^. 

If IndA = 0, then A is nonsingular, and A^ = A^^. 

Remark 2.1. Since the equation 3) of the equation 1) can be replaced by 
follows 

la) XA'^^+i A^ 

The Drazin inverse can be represented explicitly by the Jordan canonical form 
as follows. 

Theorem 2.1. 0/ // A e C"''" with IndA = k and 

where C is nonsingular and rankC = rankA*"', and N is nilpotent of order k, 
then 

We use the following theorem about the limit representation of the Drazin 
inverse. 

Theorem 2.2. 0/ // A G C"^", then 

A^ = lim (AI„ + A''^+l)"'A^ 

where k — IndA, A G M+, and M-|- is a set of the real positive numbers. 
The following theorem can be obtained by analogy to Theorem 12.21 
Theorem 2.3. If A e C"^", then 

A^ = limA'^-(AI„ + A'^-+i)"\ 
where k — IndA, A G IR.-|_. and M-|_ is a set of the real positive numbers. 
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Denote by a'-'^-' and a^^'"'' the jth column and the zth row of A*^ respectively. 



Lemma 2.1. (IMJ, Lemma 3.1) If A e C"^" with IndA = k, then for all 
i,j = l,n 



rank A 



T (a'f ) 



< rankA'=+^ 



Using Theorem 12.21 and Lemma 12.11 we obtained in [22] the following determi- 
nantal representations of the Drazin and group inverses and the identity A^A 
on i?(A'=). 

Theorem 2.4. ([2^], Theorem 3.3) If IndA = k and rankA'^+i = rankA*^' = 
r < n for A G C"^", then the Drazin inverse A^ ~ i^Fj) ^ C"^" possess the 
following determinantal representations: 



/3GJ.,„{j} 



(k) 



(5) 



for all i, j = 1, n. 



Corollary 2.1. ([22], Corollary 3.1) If IndA = 1 and rankA^ = rankA = 

r < n for A G C"^", then the group inverse A^ = (afj) G ([^"X" possess the 
following determinantal representation: 

E (A^.(a.,))^ 

„ /3GJ.,„{»} 



E (A2)^ 



for all i,j = 1, n. 

Corollary 2.2. ([22], Corollary 3.2) If IndA = k anrf rank A''+^ = rankA*" = 
r < n for A G C"^", then the matrix A^ A = {pij) G C"^" possess the following 
determinantal representation 



E (A^+^ (. 



.(fe+i) 



Pij 



)) 



(6) 



for all i,j ~ 1, n. 



Using Theorem 12.31 we can obtain another determinantal representation of 
the Drazin inverse. At first we consider the following auxiliary lemma similar 
to Lemma [2. II 
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Lemma 2.2. If A G C"^" with IndA = k, then for all i,j = l,n 
rankA*=+i (af^^ < rankA'^+i. 

Proof. The matrix A*^^^ (^j'-O represent as follows 



/ ^ (fc) 



(fc) 
Si 



s=l 



(fc) 



Let P;i (— a; j) e C"^", {I ^ i), be a matrix with —aij in the {I, i) entry, 1 in all 
diagonal entries, and in others. It is a matrix of an elementary transformation. 
It follows that 



(fe) 



'(fe) 



(fe) 



'(fe) 



(fe) \ 



(fe) 



The obtained above matrix has the following factorization. 

,(fe)\ 



(fe) 



(fe) 





n 


/an .. 





.. 


1 


\a„i . . 






,(fe) 



(fe) 



E O-nsO: 



(fe) 
sn 



Denote the first matrix by 



A := 



/ an 


\ ani 



ain\ 


^nn / 



1 



/ (fe 



(fe) Ak) 



12 

(fe) (fe) 

'^21 '^22 



. (fe) (fe) 



"in 
(fe) 
^2n 





jth 



din \ 



^nn j 



ith. 
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The matrix A is obtained from A by replacing all entries of the ith. row and 
the jth column with zeroes except for 1 in the entry. Elementary transfor- 



mations of a matrix do not change its rank. It follows that rank A*^"*"^ (^j*^') 
min rank A*^ , rank A >. Since rank A > rank A*"' the proof is completed. 



< 



Theorem 2.5. IflndA = k and rank A''+^ = rankA*-' = r < n for K G C"^", 
then the Drazin inverse — (af^) G possess the following determinantal 

representations: 



E 



A*^ 



E I (*"+'); 



(7) 



for all i,j = 1, n. 



Proof. If A G K+, then rank (AI + A'^+^) ~ n. Hence, there exists the inverse 
matrix 



AI + A*^ 



( R\\ i?21 



det (AI + A'=+i) 



In J^ln 



Rn2 
Rn n j 



where Rij is a cofactor in AI + A'''+^ for all i,j = By Theorem 12.31 

A-" = liniA*^ (AI„ + A'^+i)"\ so that 



K° = lim 



A^odet(AI + A'^'+i) 



E" „(*:) n 
Ons -Kli 



Z^s=l "is -'^"s 



lim 



det(AI+A'' 



dct(AI+A'= + i)^ 
det(AI+A'=+i) 



dot(AI+A'=+i)^ (a(^)) 
det(AI+A>fT) 



dct(AI+A'= + i)^ ) 
dct(AI+A'=+i) 



(8) 



Similar to the characteristic polynomial, we have 

det (AI + A'^^+i) = A" + diA"-i + dzA""^ + . . . + d„, 

where = E |(-^'^^^)"| ^ sum of the principal minors of A'^'"''^ of order 

s, for all s = 1, and dn — dctA'''+^. Since rankA'^+^ — r, then dn = 

dn-i = . . . = dr+i = and 



det (AI + A'^+i) = A" + diX"-^ + daA"-^ + . . . + drV' 



(9) 
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Similarly we have for all i,j = l,n 

det (AI + A'^'+i)^. (a['-)) = l['''>X"-^ + l^'^h"-^ + . . 



where for all s = 1, n — 1, 



(k) 



"6-f=,„{j} 



and l^n''' = det A 



By Lemma [2.2[ rankA*^^^ (^^'^O ~ ^' ^^^^ if s > r, then for all a £ 
Is,n{'i-} and for all i, j = 1, n, 

|(ArK')):h»- 

Therefore if r + 1 < s < n, then for all i, j — \,n, 

^ |(A,^+i(af)))].0, 

and d'^' = det A)'+^ {sl^^^^ = 0. Finally we obtain 

det (AI + A'^+i)^. (a[")) = l^^'^X'-^ + 4'^'^A"-^ + . . . + ^fe)A"-^ (10) 

By replacing the denominators and the nominators of the fractions in the 
entries of the matrix (|H]) with the expressions © and dTU]) respectively, finally 
we obtain 



A^ = lim 

A-i-O 



A"+diA"-i + ...+drA"- 

i<"^A"-i + ...+i("i)A"- 
■ A'>+(ii A"-i + ...+drA"- 



;(1")a"-i + ...+;(1")a" 
A'»+£iiA"-i + ...+drA" 



i<"">A"-i + ...+i^"">A" 
A"+diA"-i + ...+(i,-A" 



where for all i, j = 1, n, 

0^ 



E 

Q6-fr.„{i} 



, fc+1 



This completes the proof. 
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Using Theorem [23] we evidently can obtain another determinantal represen- 
tation of the group inverse and the foUowing determinantal representation of 
the identity AA^ on i?(A'=) 

Corollary 2.3. If IndA = 1 and rankA^ = rank A = r < n for A e C"""", 
then the group inverse A^ = (afj) G possess the following determinantal 

representations: 



E i(A2):i 



(11) 



for all i,j = 1, n. 



Corollary 2.4. If IndA — k and rankA'^^^ = rank A''' = r < n for A e 
C"^", then the matrix A A^ = {q^j) £ C"^" possess the following determinantal 
representation 



q^J = 



(12) 



for all i,j = 1, n. 



3. Cramer's rule of the Drazin inverse solutions of some matrix equa- 
tions 



Consider a matrix equation 



AX ^ B, 



(13) 



where A G C"^" with IndA = k, B (E C"^" are given and X e C"^" is 
unknown. 

Theorem 3.1. (fl^], Theorem 1) If the range space i?(B) C i?(A'=),, then the 
matrix equation Iil3\) with constrain i?(X) C R{A'') has a unique solution 

X = A^B. 

We denote A'^B =: B = (%) G C"^". 

Theorem 3.2. //rankA''+^ = rankA'' = r < n for A E C"^", then for 
Drazin inverse solution X — A-^B = (xij) G C"x™ qJ {jgj) we have for all 



i = '^,n, j = l,m, 



E 

PeJ,■.,^{i} 



E (A'^+i) i 



(14) 
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Proof. By Theorem 12.41 we can represent the matrix by Therefore, 
we obtain for all i = l,n, j — 1, m, 



E ELr (A^rfa'.?)) 



E i^'+')^B 



Since X] ^''■s^ 



b then it follows (fT4l). 



Corollary 3.1. (^/i?^/, Theorem 4-2.) IflndA ~ k and rank A''+^ = rankA'^ = 
r < n for A £ C"^", and y ~ (yi, . . . ,yn)"^ G C", t/ien /or Drazin inverse 
solution X = A^y =: . . . , a;„)"^ £ C" o/ f/ie system of linear equations 



we have for all j ~ \,n, 



A-x = y, 



E__ {^Tnl 



(A'^+i)^ 



where f = A'^y. 

Consider a matrix equation 

lere A 
unknown. 



XA = B, 



(15) 

are given and X G C"^™ is 



Theorem 3.3. (^], Theorem 2) If the null space 7V(B) D iV(A''), then the 
matrix equation !il5\) with constrain -/V(X) D N{A^) has a unique solution 



X ^ BA^ 



We denote BA*^ =: B = (%) e 
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Theorem 3.4. //rankA'^+^ = rankA'^ = r < m for A e C™^™, then for 
Drazin inverse solution X = BA"^ = i^ij) £ C"^™ of HIS]) , we have for all 



i^^,n, j = l,m, 



k+l 



E l(A^-+^)gl 

QG/,. m 



(16) 



Proof. By Theorem 12.51 we can represent the matrix A^ by ([T]). Therefore, 
for all i = 1, n, J = 1, m, we obtain 



E l(A'=+^)gl 



E l(A'^-+^)2l 



Since for alH = 1, 



(fc) 



ism 



then it follows ([T| 

Consider a matrix equation 



AXB = D, 



(17) 



where A e C"^" with IndA = fci, B G 
are given, and X g C"^™is unknown. 



^"^ with IndB = fc2 and D £ C"^'" 

Theorem 3.5. ((JJ], Theorem 3) If R{Ti) C R{A^^) and iV(D) D N{B'"'), 
k = max{ki, ^2}, then the matrix equation |_?9| ) ivit/i constrain i?(X) C i?(A'^) 
and -/V(X) D A'^(B*'') /las a unique solution 

X = A^DB^. 

We denote A'^^DB'^-^ =: D = (d^) G C"^™. 

Theorem 3.6. // rankA'^i+i = rankA^^^ = ri < n /or A G C"""", and 
rankB'^^"*'"'^ — rankB*^'^ — r2 < m for B G C™^"\, then for the Drazin inverse 
solution X = A^DB-" =: {x^j) G C"^"' 0/ {I^j we have 



E 



E (A'^^+i)^ E |(B'=^+i):|' 



(18) 
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or 



E 

ie/r2,™{j} 



B 



k2+l 



(19) 



where 



(a.) 



aelr2.m{j} 



-I T 



^(d„.) 



(20) 



are the column-vector and the row-vector, dt, and d,j are respectively the ith 
row and the jth column of D for all i = \,n, j = 1, m. 

Proof. By Theorems HH] and [13] the Drazin inverses = (a^) e 
and = (6g) G 



respectively, 



mxm possess the following determinantal representations, 



a. 



E (A^^+^) ^ 



E B^^^+i(bf^))2 



a6/r2.m{i} 



E l(B'=^+i)g| 



(21) 



Then an entry of the Drazin inverse solution X = A^DB-^ =: {xij) G C"^"^ is 

m / n \ 

Xii=Y.[Y.a^d,Ah%. (22) 



s=l \t=l 



Denote by d.^ the sth column of A*D D = (dy ) G C"^™ for all s = l,m. 
It follows from ^ a^fdts = d. s that 



E 

/3eJri,„{i} 



i=l 



t=l 



E 



E E A'^^.^^ (a^^)) 
/3eJ.i.„{i}t=i ^ ^ 



dts E AY^fd.) 



E (A'^-^+i) 



E (A'=^+i) ^ 



(23) 
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Substituting and (gTI) in (HH), wc obtain 



E 



^ E (A'=^+^) 



Suppose e^. and e. ^ are respectively the unit row-vector and the unit column- 
vector whose components are 0, except the sth components, which are 1. Since 

n rn m 



1=1 



s=l 



then we have 



m m n 



EEE E A'^+'iB.,)', 



=lt=l/=l/3GJ^i,„{!:} 



E B^^+i(e02 



Qe-fr2,m{j} 



E (A''-^+^);^ E l(B'=^+^)g| 



EE E 

t=i i=i ^eJri, 



,fcl+l 



(e.O 



ae/r2,7Ji{i} 



E (A'=^+^)^ E l(B''^^+^)2l 



Denote by 



4 



(24) 



the tth component of a row- vector = (d^, df^) for aU i = 1, m. Substi- 
tuting it in ([24| . we obtain 



Erf^ E B^^+i(eOS 



*=1 ae/,.2,„{i} 



E (A'^^+^)^ E l(B''^^+^) 



Since ^ ^^(64. — df-, then it fohows ((T 
t=i 

If we denote by 



dfr--Y.dit E Ib^'K)" = E |Br'(dOS 
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the Ith. component of a column- vector = (dy, for all I = l,n and 

substitute it in (p4|) . we obtain 



Since V e /d^ = d^,, then it follows (HH 



4. Applications of the determinantal representations of the Drazin 
inverse to some differential matrix equations 



Consider the matrix differential equation 

X' + AX = B 



(25) 



where A G C"^", B G C"^" are given, X e C"^" is unknown. It's well-known 
that the general solution of (|25|) is found to be 



X(i) = exp 



-At 



exp^* dt \ B 



If A is invertible, then 



exp^* dt = exp^' +G 



-1 „„„At 



where G is an arbitrary n x n matrix. If A is singular, then the following 
theorem gives an answer. 



Theorem 4.1. ([26], Theorem 1) If A has index k, then 



exp^* dt ^ A^ exp^' +(I - AA^)i 



. D At 



A* 



k~l 



fc! 



Using Theorem 14.11 and the power series expansion of exp we get an explicit 
form for a general solution of (1^51) 



X(t) = 

{A^ + (I-AA^)t(l 



2 i + 31 t 



.(_l)fc-lAl^t/c-l^ ^q|b^ 

If we put G = 0, then we obtain the following partial solution of ([25]) . 



X(i) = A^B + (B - A^AB)t - i(AB - A^A^B)t^ 



i-i)" 



(A*-iB- A^A'^-B)i*. 
Denote A'B =: B^') = (b[^h e C"^" for all / ^ T^. 



(26) 
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Theorem 4.2. The partial solution i26]) . X(t) — {xij), possess the following 
determinantal representation, 



E |(A'=t^(b(;+^))) 

1 I 5(1) _ ' ^ ^ / / 



3e Jr. n 



+ 



Mr 



r(fc-l) geJr,,i{.} 



(A^ + i(b<f)))|| 



E |(A'^'+^)^| 



(27) 



for all i,j = 1, n. 

Proof. Using the determinantal representation of the identity ^ we obtain 
the following determinantal representation of the matrix A^A"'B :— (yij), 



s=i t=i /3eJr.„{«} 



F — ■ 



•6* 



E 

/3ej,,„{i} 



l3£Jr.n 



E (A'^+i) ^ 

/3GJr. 



for all z, J = 1, n and m = 1, fc. From this and the determinantal representation 
of the Drazin inverse solution (fT4|) and the identity ([6|) it follows ([27|) . 



Corollary 4.1. If IndA = 1, then the partial solution of ^¥5\ 

X(t) = (xy ) = A^B + (B - Af AB)i, 
possess the following determinantal representation 



E (A2) ^ 



— + h - - 

V 



for all i,j = 1, n. 

Consider the matrix differential equation 

X' + XA = B 



E (A2) ^ 



i. (28) 



(29) 
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where A G C"^", B G C"^" are given, X G C"^" is unknown. The general 
solution of (1291) is found to be 



X{t) =Bexp-^* (^J exp^* 



dt 



If A is singular, then an exphcit form for a general solution of (|29|) is 

m = 

B {a^ + (I - AA^)t (l - #t + + ...(_i)fc-i Al^t^c-i^ ^ q| 

If we put G = 0, then we obtain the following partial solution of ([29l) . 
X(i) = BA^ + (B - BAA^)t - i(BA - BA^A^ji^ + ... 



i-jr 



'{r-CBA"-^ - BA'=A^)i'= 



(30) 



Denote BA' B^) = {b[^^) G C"^" for all I = 1,2k. Using the determinantal 
representation of the Drazin inverse solution (fT6|) . the group inverse ([TT|) and 
the identity (fT2|) we evidently obtain the following theorem. 

Theorem 4.3. The partial solution i30\} . X(t) — {xij), possess the following 
determinantal representation, 



^1.3 



2 I E l{A'=+i)3t 

E KAj+Vbf^')) 
(-1)" I _ °g^^.>»{^>' ^ ^ ^ ■ 

fe! I 



E l(A'=+i)g| 



t2 + 



E l(A'=+i); 



/or a// «, J = 1, n. 

Corollary 4.2. If IndA — 1, then the partial solution of Ii29\) . 

X(t) = (x^j) = BAf + (B - BAA3)t, 
possess the following determinantal representation 



E l(A2)g| 



E l(A2)Sl 



/or all i,j — 1, n. 
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5. Examples 

In this section, we give examples to illustrate our results. 

5.1. Example 1 

Let us consider the matrix equation 



AXB = D, (31) 



where 







i 









i 




i 





B = 

We shall find the Drazin inverse solution of ([31]) by ([T8|) . We obtain 

4 ON 
= I 2 - 2i 

-2-2i Oy 



—I 

B'=\ 1 















4-4i 










\-4 - Ai 








i 


3- A 






-1 


1 + 3i . 






i-i 


3 + 1 







Since rank A = 2 and rankA^ = rankA'^ = 1, then ki — IndA = 2 and ri = 1. 
Since rankB — rankB-^ = 2, then fc2 = IndB — 1 and r2 — 2. Then we have 



and ^ (A3)'^ 

/3eJi.3 

E |(B^ 

"€^2, 3 



4-4i 
-4-4i> 




+ (-9 - 9i) + (9 - 9i) = -18i. 
By ([20]) . we can get 



= I -12i I , = I 12i ], d^^ 
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Since A\ {d^^) 



'12 - 12i 0^ 

— 12i I , then finally we obtain 

-12 0; 



Xll 



E A3,(dB)^ 

/9eJi.3{i} 



12-12i l + i 



Similarly, 



Xl2 = 



-12 + 12i 



8 ■ (-18i) 



12 



Xl3 



(-18i) 18' 



-12i 



3^21 



8-(-18i) 12 



2^22 



12? 



-12-12i 1-i 



(-18*) 



12' 8 • (-18i) 12 



12 



-12 



a;3i = 



Then 



8 • (-18«) 



12' 



8 • (-18i) 



-12 + 12i 
8- (-18i) 





-l-i 


2 


12 


12 




1 


1 


l — l 


12 


12 


12 


i 


i 


-l-i 


12 


12 


12 


(ED. 







-l-i 

12 ■ 



5.^. Example 2 

Let us consider the differential matrix equation 



where 



X' + AX = B, 



B 





(32) 



Since rank A = rankA-^ = 2, then k = IndA = 1 and r — 2. The matrix A is 
the group inverse. We shall find the partial solution of by ([25]). We have 

-i i /2-i 2i 0^ 

A^=| 1 -1 l + 3i , B(i) = AB= l + 2i -2 oj, 



'2-2i 2 + ii 0> 
§(2) ^ A^B = I 2 + 2i -i + 2i 
,l + 5i -2 0> 
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and 



as J2, 

det 



.2^ /3 



det 



-1 l + 



1 -IJ 3 + i 

0+(-9-9i) + (9-9?:) = -18i. 



det 



— i 3 — i 
+ i 3 + i 



Since (A2) ^ (b^J^) 



2-i 



then finally we obtain 





3 - 


M 






1 + 


3i and 




i 


34 


V 




2 


- 2i 


i 3 


— i 




+ 2i 


-1 H 


-'Si 




+ 5i 


3-i 3 


+ i 



Xii 



E |(A^(b<;)))|| 



E |(A^(b(?'))g| 



E |(A^)^| 



fie J2,3 



3-3i 
-18i 



/ — 1+i 
-18t I 6 • 



Similarly, 



X22 



-3 + 3i 
-ISi 



-3-3i 
-ISz 



9 + 9A 1 + i ^ , 

« — I < = + ^^t, a:i3 = + (1 



2:21 = 



-18i 

3 + 3i 
487 



0)i = i, 



-9 + 9i' 



-18\ -1 + i 

I t = , 

-18iy 6 ' 

1-i 1+i 



g +^^' a;23 =0 + (l 



0)t = t, 



X31 



X32 



9 + 3i 
-18i 



Then 



X = - 
6 



-12i 



-18i 
-18 



-18A 2 
1 ]t= -, 

-18i 3' 



-18i 



-l + 3i , , 

t= , X33 -0+(0-0)t = 0. 



1 + i -l-i + (3 + 3i)i t 
-1 + i l-i + (3 + 3i)t t 
4 -l + 3i 



is the partial solution of 
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